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In this paper we present detailed Mossbauer spectroscopy study of structural and magnetic properties 
of the undoped parent compound CaFe2As2 single crystal. By fitting the temperature dependence of 
the hyperfine magnetic field we show that the magneto-structural phase transition is clearly first-order in 
nature and we also deduced the compressibility of our sample to be 1.67 x 10~ 2 GPa -1 . Within the 
Landau's theory of phase transition, we further argue that the observed phase transition may stem from 
the strong magneto-structural coupling effect. Temperature dependence of the Lamb-Mossbauer factor 
show that the paramagnetic phase and the antiferromagnetic phase exhibit similar lattice dynamics in high 
frequency modes with very close Debye temperatures, 0£> ~270K. 

PACS numbers: 76.80.+y, 74.10.+V 



INTRODUCTION 

The discovery of superconductivity (SC) with critical 
temperature (T c ) up to 55 K in iron-arsenide systems 0]- 
3] has triggered enormous interest in iron-based supercon- 
ducting materials. Generally, the parent compound un- 
dergoes a structural transition and exhibits antiferromag- 
netic (AFM) order below room temperature. Doping with 
electrons or holes into the parent compound will suppress 
the structural and magnetic transitions and induce SC 0]. 
Nowadays it is widely believed that SC is directly cou- 
pled to the magnetism in iron-based superconductors JHQ]. 
There is also clear evidence about the strong spin-lattice 
coupling 171 [8fl. Meanwhile, though the electron-phonon 
(e-p) coupling alone can not explain the high critical tem- 
perature, it is proved that the e-p couplingeffect definitely 
plays some role in the SC mechanism ]9|, [lOj], probably 



through the spin-channel I1UI12H . Hence a better under 



standing of the magneto-structural coupling effect and the 
lattice dynamics in the iron- arsenide parent compounds is 
important in unraveling the SC pairing mechanism of the 
iron-based superconductors. 

57 Fe Mossbauer spectroscopy (MS) is an excellent probe 
for both structural and magnetic local properties in iron- 
containing compounds. In iron-based superconductors 
the iron element serves as the primary constituent, which 
makes MS a useful analytical tool 1131] . especially for in- 
vestigating the lattice dynamics in these systems. Indeed, 
MS has been widely used to study the iron-based supercon- 
ductors, e.g. ReFeAsO (Re=rare earth elements) Ill4 Il5ll . 
AFe 2 As 2 (A=Ba, Ca, Eu) system and LiFeAs 

lfl9h . However, properties regarding magneto-structural 
coupling effect and lattice dynamics are hardly discussed 
in most of these works. In the present work, single crys- 
tals of CaFe 2 As 2 were synthesized and studied in detail by 
MS in the temperature range of 18 K to 290 K. It is found 
that upon cooling through 170 K, CaFe 2 As 2 undergoes a 
first-order magneto-structural phase transition and changes 



from a tetragonal ThCr 2 Si 2 -type paramagnetic (PM) phase 
to an orthorhombic AFM phase B20U21I1 . Our observations 
suggest that the lattice dynamics of both the tetragonal and 
orthorbombic phases are very similar in the high-frequency 
modes, which could be denoted by the almost same Debye 
temperatures of the two phases. 



EXPERIMENTS 

Single crystals of CaFe 2 As 2 were grown by the FeAs 
self-flux method described elsewhere lEfjIl . The resis- 
tivity and hall effect analysis were reported earlier lEfjIl . 
The phase purity of the crystalline sample was confirmed 
by single crystal X-ray diffraction measurements using a 
Philips X'pert diffractometer with Cu K a radiation. DC 
magnetization measurement was carried out through a 
commercial (Quantum Design) superconducting quantum 
interference device (SQUID) magnetometer. The transmis- 
sion Mossbauer spectra were recorded using a conventional 
constant acceleration spectrometer with a 7-ray source of 
25mCi 57 Co in palladium matrix moving at room tem- 
perature. The absorber was kept static in a temperature- 
controllable cryostat. All isomer shifts are quoted relative 
to a-Fe at room temperature. The Mossbauer spectra are 
fitted with MOSSWINN 3.0l programme JH]. 



RESULTS AND DISCUSSION 

Fig- CD presents the single crystal X-ray diffraction pat- 
tern of the CaFe 2 As 2 compound. The inset shows the tem- 
perature dependence of the magnetic susceptibility mea- 
sured under H = 5 T. As can be seen, only (001) diffraction 
peaks are observed, indicating that the crystallographic c- 
axis is perpendicular to the plane of the plate -like single 
crystal. The susceptibility curve shows an anomaly around 
165 K corresponding to the magneto-structural phase tran- 



sition, which is similar to previous reports [20]. All these 
results indicate that the sample has a good quality and 
hence a good starting point for the following MS measure- 
ments. 
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FIG. 1. Single crystal X-ray diffraction pattern of CaFe2As2. The 
inset shows the temperature dependence of the magnetic suscep- 
tibility measured under H = 5 T. 

The 57 Fe Mossbauer spectra of the single crystal 
CaFe 2 As 2 recorded in the temperature range of 18 K-290 K 
are shown in Fig. |2j Above 170 K, the spectra are 
quadrupole split paramagnetic lines and can be well fit- 
ted with only one asymmetric doublet. This means that 
the local environment of the Fe ion is unique, indicating 
no impurity phases exist in the Mossbauer absorber. The 
spectra recorded at 160 K and 165 K were fitted with a sex- 
tet corresponding to the AFM phase and an asymmetric 
doublet corresponding to the residual paramagnetic phase. 
The coexistence of both phases is characteristic for first- 
order phase transitions. At temperatures below 155 K, the 
spectra can be fitted with the superposition of two Zee- 
man patterns. Two subspectra are needed to get a better 
fit to the observed spectra, which is similar to previous re- 
ports [23]. This may be understood as both commensu- 
rate (main component) and incommensurate (minor com- 
ponent) spin density waves (SDW) exist in the sample. The 
incommensurate SDW has also been observed in Ref i23"ll . 
where three sextets are used to fit the Mossbauer spec- 
trum. Yet the incommensurability in 122-type iron-based 
parent compounds seems to be controversial. For exam- 
ple, for BaFe 2 As 2 and EuFe 2 As 2 compounds one compo- 
nent is enough to get a good fit to the spectrum and the 
shape of the SDW is quasirectangular [17] at low tempera- 
tures. While the shape of the SDW does not develop a fully 
quasirectangular shape Ill7ll even at liquid helium temper- 
ature, which may be due to the incommensurate SDW. 
So, understanding the incommensurate SDW certainly de- 
serves much more careful experimental examinations. 

The temperature dependence of hyperfine parameters ex- 
tracted by fitting the measured Mossbauer spectra at dif- 
ferent temperatures are plotted in Fig [3] It is well known 
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FIG. 2. 57 Fe Mossbauer spectra of the single crystal CaFe 2 As 2 
obtained at indicated temperatures. The spectra taken between 
18K-155K were fitted with two Zeeman patterns as discussed 
in the text. The spectra recorded at 160 K and 165 K were fitted 
with a sextet and an asymmetric doublet. The spectra taken above 
170 K were fitted with only one asymmetric doublet. 



that the total isomer shift (5) is the sum of chemical shift 
(5c) and thermal shift (5 S od), i-e. 5 = 5 C + 5sod- As 
shown in Fig [3] (a), the isomer shift has a typical value 
of 0.42 mm/s for the 122-type iron-based superconduc- 
tors at room temperature 11311 . Upon cooling, 5 increases 
gradually before 170 K, the characteristic of 5 sod, then 
shows a sudden increase in a narrow temperature range, 
which could be mainly ascribed to the contribution of 5q. 
In Mossbauer experiments, 5c can be expressed as 112411 
5c = ct(p a Fe — Pa), where p a Fe and Pa denote the elec- 
tron density at the iron nucleus for a-Fe and the absorber, 
respectively, a is a calibration constant, which is positive 
for 57 Fe nucleus. Thus, a positive isomer shift indicates a 
smaller electron density at the iron nucleus in the sample 
than in a-Fe. Therefore, the sharp increase in 5 at ~165 K 
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corresponds to a decrease in the electron density at the iron 
nucleus, which may arises from the changes of the unit vol- 
ume at the transition temperature. The structural transition 
results in larger Fe-As distances in the orthorhombic phase, 
which reduces the interaction between Fe-3d and As-4p 
electrons and leads to more localized Fe-3d electrons. This 
will enhance the screening effect of the Fe-3d electrons to s 
electrons near the nuclear, causing smaller s-electron den- 
sities at the iron nucleus thus a bigger isomer shift. 
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FIG. 3. Temperature dependence of the Mossbauer hyperfine pa- 
rameters extracted by fitting the spectra taken at different tem- 
peratures, (a) Isomer Shift 5, (b) Line Width T, (c) Quadrupole 
splitting | A.Eq \. 

The spectrum line width, T, and quadrupole splitting, 
| AEq |, are shown in Fig[3](b) and (c), respectively. They 
both show an anomaly around 165 K, which is ascribed to 
the magneto-structural phase transition. One should no- 
tice that the line broadening actually appears at temper- 
atures above the phase transition. Similar line broaden- 
ing effect has been observed in (Bai_ x K x )Fe2As2 systems 
|25J. The most probable reason for the line broadening 
effect can be considered as the onset of magnetic order 
and/or short range magnetic fluctuations that presents in the 
CaFe 2 As 2 compound [8]. These observations highly sup- 
pose the magnetic origin of the magneto-structural phase 
transition in CaFe 2 As 2 parent compound. Upon cooling, 



some Fe spins happen to form clusters with short-ranged 
magnetic order in the nonmagnetic matrix by spin fluctu- 
ations of the Fe-3d charges. These spin clusters broaden 
the line width instead of detectable hyperfine field. With 
further cooling, the phase transition happens when the spin 
clusters reaches its critical point and causes structural dis- 
tortion simultaneously to optimize the energy of the sys- 
tem due to spin-lattice coupling effect. Therefore, it is sug- 
gested that the magnetic fluctuations and spin-lattice cou- 
pling effect may be responsible for the magneto-structural 
phase transition in CaFe 2 As 2 compound. 



l 

13 

Sh 0.8 

a 

EH 

<3 0.6 
& 

X 

> 04 
td 

1 0-2 









B u i0) = 9.73 T 
hf 




T Q = 170 K 






Hi • 



0.2 0.4 0.6 0.8 

Relative Temperature 



1.2 



FIG. 4. The reduced hyperfine magnetic field vs reduced temper- 
ature for the single crystal CaFe 2 As 2 sample. The solid line is 
theoretical fit to the experimental data obtained on heating, as- 
suming Tb=170 K (see text). 

To examine the magnetic properties and get insight into 
the AFM phase transition, the temperature dependence of 
the hyperfine magnetic field is investigated using a mean 
field model after Bean and Rodbell 12611 . which, while de- 
scribing the order parameter variation correctly, does not 
yield information about the magentic fluctuations in side 
of the hysteretic region, observed here. The reduced hy- 
perfine magnetic field vs reduced temperature is plotted in 
Fig H] and the solid line is the calculated theoretical curve. 
In the Bean-Rodbell model, the temperature dependence 
of the sublattice magnetization in a AFM material is ex- 
pressed as [26] T/T = [a I tanh -1 cr)(l + r/a 2 /3) in the 
low spin state approximation, r] = ^Nk B KT Q f3 2 , where 
a is the reduced sublattice magnetization, r] is a fitting pa- 
rameter (rj < 1 for a second order transition, r] > 1 for 
a first order transition and for ij = the equation reduces 
to a Brillouin function), /3 is the slope of the dependence 
of transition temperature on volume, T would be the tran- 
sition temperature if the lattice were not compressible, K 
is the compressibility, T is the temperature and k B is the 
Boltzmann constant. As can be seen, Fig |4] shows an ap- 
proximate agreement between the theory and experimental 
data. The observed hyperfine field decreases more faster 
with increasing temperature than the theory predicts in the 
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low temperature range, which could be explained by the 
low-lying spin excitations as evidenced by inelastic neu- 
tron scattering [8]. The fitted values of rj and B h f(0) with 
fixed T =170K are found to be 1.35 and 9.73 T, respec- 
tively. The value of r] indicates that the nature of the AFM 
transition is of first order. This is in agreement with previ- 
ously reported thermodynamic, transport and microscopic 
data Il20n . And the fitted zero point hyperfine field is very 
close to the value (B h f^lOT) reported by Kumar et al. 
ill at 4.2 K. 

If we assume the Bean-Rodbell model to be valid, we 
may determine the compressibility of the CaFe 2 As 2 com- 
pound. According to the Bean-Rodbell model, the com- 
pressibility is given by K = 6(Av/v ) 2 / (rjNkBT a 4 ), 
where Av/vq is the volume change at the transition. If 
we take from experimental values, the volume change 
Av/v = —0.93% 12711 (The negative sign of Av/v 
means a compression of the lattice during the AFM/PM 
transition, which is consistent with the Bethe-Slater 
curve.), T = 170 K, a = 0.77 (taken at T=165 K), rj = 
1.35 and, assuming N = 2.79 x 10 28 m~ 3 (corresponding 
a density of 7 g /cm 3 ), one gets K = 1.67 x 10~ 2 GPa~ x . 
The results is very close to the value of first principles 
calculations on CaFe 2 As 2 (K = 1.65 x 10~ 2 GPq- 1 ) 
1 2811 . while a little larger than that of the reported experi- 
mental values, 1.04 x 10~ 2 GPoT 1 for the LaFePO El 
and 0.98 x 10~ 2 GPa' 1 fortheNd(0 . 8 8Fo.i2)FeAs jSj. 
Thus a direct measurement of the compressibility would be 
a useful check on this theory. 

As discussed above, the magneto-structural phase tran- 
sition in CaFe 2 As 2 is first-order in nature. In order to 
obtain a theoretical interpretation of our experiments, we 
adopt the phenomenological Landau theory of phase tran- 
sitions to the system under study. In Landau's theory the 
free-energy density of the system is expanded in terms of 
two order parameters of a (sublattice magnetization) and e 
(strain) as JH 

F(a, e) = -a(T)a 2 + \ba 4 + -ca e + -Be 2 + Xe ■ 0) 
2 4 6 2 

Where B is the elastic modulus and A is the magneto- 
elastic coupling coefficient. Minimizing the total energy 
with respect to the strain yields e = —\a 2 /B. Substitut- 
ing e with a in Eq. (TJ), one obtains a renormalized 2-4-6 
Landau free energy, 

F(a) = -a(T)a 2 + (-b - — )cr 4 + -ca e . (2) 

The magneto-elastic coupling is presented by the fourth or- 
der term. To estimate the coefficient of this term, we turn to 
the Bean-Rodbell model, from which the magnetic part of 
the free energy without external field is also a 2-4-6 func- 
tion of a, as §32 



Comparing Eq. (O and (0, clearly, the sign of (6/4 — 
X 2 /2B) is determined by (T — r]T ), which is negative 
when r] = 1.35 around the transition temperature. A neg- 
ative fourth order term in Landau free energy creates an 
energy barrier in the free energy landscape, this leads to 
a first-order transition. Obviously, a large magneto-elastic 
coupling coefficient A results in a large negative fourth or- 
der term. The more negative is this term, the larger is the 
transition barrier, and the more distinct are the first-order 
transition characters. On the other hand, a larger A leads 
to a larger spontaneous lattice distortion upon the mag- 
netic transition. Therefore, the strength of the magneto- 
elastic coupling A can be represented by the magnitude of 
the spontaneous lattice distortion, like rj, which involves 
the parameters K and (3 that are related to the volume 
change, rj also controls the order of the magnetic phase 
transition. According to Eq. (0, one always gets a nega- 
tive fourth order term when rj > 1, indicating a first-order 
magnetic phase transition at T due to the energy barrier 
from magneto-elastic coupling. Thus, the parameter r\ pre- 
dicts the strength of the coupling. 

To investigate the lattice dynamics of the Fe atoms, the 
temperature dependence of the Lamb-Mossbauer factor or 
recoil-free fraction, /, is discussed in what follows. In De- 
bye model of the lattice dynamics, / can be expressed as 



f(T) 



T,2 r<3> D /T _ 



- bt -R fl I I T -.2 K- 
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(4) 



where k B is Boltzmann's constant, ©£> is the Debye tem- 
perature, and E R is the recoil energy of an free emitting 
nucleus. In high temperature (T > Q D /2) approximation, 
we may rewrite Eq. © as the following formula 12411 : 

-Hf(T)) = + eT + 5T 2 + •••). (5) 

Here e and 5 are termed as the anharmonic coefficient and 
they are both very small values. 




G{a) = \{T~ T )a 2 + ±(T - r?T )a 4 + ±ca 6 .(3) 
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FIG. 5. Logarithm of the relative / factor versus temperature. 

The temperature dependence of —ln(f(T) //(18 K)) is 
shown in Fig [5j It is evident that upon cooling, the curve 
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drops steeply through 170 K to 150K due to the occur- 
rence of the phase transition. This sudden increase in f(T) 
is likely due to the thickness effect, i.e. the increase of 
spectral area when overlapping peaks separate from each 
other, as it usually happens when a magnetic phase transi- 
tion takes place [22]. The temperature dependence of f(T) 
beyond this temperature region is quite linear, i.e. Fe atoms 
experience harmonic vibrations in both tetragonal phase 
and orthorhombic phase. Since the anharmonic effect char- 
acterized by e and 5 is negligible, the f(T) data can be sim- 
ply fitted with -Zn(/(T)// (18 K)) = 6E R T/k B Q 2 D + 
ln(f (18 K)), where Zn(/(18 K)) was treated as a free pa- 
rameter. The best-fit of the experimental data shown in Fig. 
[5] yields the following results: Q D = 272 K for the seg- 
ment below the phase transition and Q D = 271 K for the 
segment above the phase transition. The nearly same be- 
haviors of f(T) with temperature indicates that no obvious 
changes of the phonon modes in high frequency ranges of 
the Fe vibrational state, which is in agreement with inelas- 
tic neutron scattering results [33] on the phonon spectrum. 
Moreover, the observed harmonic vibrational behavior of 
the iron atoms indicates that the large phonon linewidths in 
CaFe 2 As 2 observed by neutron scattering 13411 may be due 
to to strong electron-phonon coupling effect. Furthermore, 
theoretical calculations of the phonon spectrum show that 
including magnetism does not improve agreement with ex- 
periments [34] for CaFe 2 As 2 . This is rather different for 
other 122-type, SrFeAsF Q and LaFeAsO [36] 1111- 
type systems, where the authors show that including mag- 
netic degree of freedom greatly improves the agreements 
with experiments. Thus more work is necessary to under- 
stand these properties. 



CONCLUDING REMARKS 

To conclude, we have studied the structural and mag- 
netic properties of CaFe 2 As 2 single crystal in detail by 
Mossbauer spectroscopy. Our results show that the struc- 
tural distortion and magnetic transition are coalesced into 
one single first-order magneto-structural phase transition 
through strong magneto-structural coupling effect. Short- 
range magnetic order and/or spin fluctuations appears well 
above the magneto-structural phase transition temperature 
and may be responsible for the phase transition. Ana- 
lyzing the temperature dependence of the hyperfine mag- 
netic field in the Bean-Rodbell model, we are able to de- 
termine the compressibility of the CaFe 2 As 2 crystal, K = 
1.67 x 1(T 2 GPa- 1 . Both the PM and AFM phases ex- 
hibit similar lattice dynamic behaviors of the Fe atoms in 
high-frequency modes and they both have a Debye temper- 
ature close to ~270 K. 
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